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Abstract 

We introduce an operational approach to the use of pulsating sources, located at 
spatial infinity, for defining a relativistic positioning and navigation system, based 
on the use of four-dimensional bases of null four-vectors, in flat spacetime. As a 
prototypical case, we show how pulsars can be used to define such a positioning 
system. The reception of the pulses for a set of different sources whose positions in 
the sky and periods are assumed to be known allows the determination of the user's 
coordinates and spacetime trajectory, in the reference frame where the sources are 
at rest. We describe our approach in flat Minkowski spacetime, and discuss the 
validity of this and other approximations we have considered. 

Key words: relativistic positioning system, astrometry and reference systems, null 
frames. 



1 Introduction 



The current position ing systems, such as GPS and GLONASS (lAshbyl . l2003l : 
Pascual-Sanchea . 120071 ) . are essentially conceived as Ptolemaic, since they are 
based on a reference frame centered in the Earth, and Newtonian, since po- 
sitioning is defined in a classical (Euclidean) space and absol ute tirne, ove r 
which relativistic (post-Newtonian) corrections are introduced (lAshbyl . 120031 ) . 
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Furthermore, these systems are effective for positioning on the Earth, but they 
can hardly be used for navigating in the space o utside the Earth, as in the Sol ar 



System and beyond. In contrast , some authors (IColll.l2006l:IColl et al.l.l2006alJb 



Rovelli I2OO2I : iBlagoievic et all . 120021 : iRuggiero fc Tartaehai 120071 : iBini et al. 



20081 ) recently proposed to use the worldlines of electromagnetic signals, emit- 
ted by objects in geodesic motion, in order to build a relativistic positioning 
system, based on the use of the so called emission coordinates. These relativis- 
tic positioning systems are also autonomous or autolocated, since any user can 
determine its own position (and spacetime trajectory) by solely elaborating 
the received signals. 

The simplest way of understanding what emission coordinates are, is to con- 
sider four emitting clocks, in motion in spacetime, broadcasting their proper 
times: the intersection of the past lightcone of an event with the worldlines 
of the emitting clocks corresponds to the proper times of emission along the 
worldlines of the emitters; these proper times are the emission coordinates 
of the given event. For example, one may think of a set of satellites orbiting 
around the Earth and equipped with onboard clocks, however, such a system 
would hardly be effective for the navigation in the Solar System ; for that pur - 
pose, a set of pulsars could rather be used fjCoU fc Tarantolal . |2004| . |2009| ). 
In fact, known pulsars emit their signals at a highly r egular rate (th i s is t he 
case, in particular, of the millisecond pulsars, see e.g. iKramer et al.l (120041 )). 
which makes them natural beacons for building a relativistic positioning sys- 
tem. What can be measured with great accuracy is the arrival time of the N-th 
pulse, so that counting these pulses can in principle allow to define something 
similar to the emission coordinates, even though in this case the actual proper 
emission time is unknown and inaccessible. 



Actually, the idea of using pulsars as st ellar beacons has been considered 
since the early years of dis covery of pulsars (ll)ownsl.ll9741) . and o t her p roposals 
are actually under study ( Sala et al. . 2004; Sheikh et al. . 20061 . 2007 ). which 
are based on the accurate measurement of the Times of Arrival (TOA) of the 
pulses or the phase differences between pulses: in general, in these approaches 
positioning is not autonomous, but can be referred to a reference frame cen- 
tered at the Solar System barycenter (SSB), because pulses arrival times or 
their phases have to be related to their expected values in the SSB. 



However, thanks to the use of emission coordinates (jCoU fc Tarantolal . 120091 ) , 
autonomous positioning by means of pulsars can be defined: in this context, 
we introduce here an operational approach to the use of periodic so urces (as- 
sumed to be at rest at spatial infinity) such as pulsars (anticipated in lTartaglia 
(120091 )). for defining an autonomous relativistic positioning and navigation 
system, based on the use of a four-dimensional basis of null four-vectors. We 
assume that a user is equipped with a receiver that can count pulses from a set 
of sources whose periods and positions in the sky are known; then, reckoning 
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the periodic electromagnetic signals coming from (at least) four sources and 
measuring the proper time intervals between successive arrivals of the signals 
allow to localize the user, within an accuracy controlled by the precision of the 
clock he is equipped with. This system can allow autopositioning with respect 
to an arbitrary event in spacetime and three directions in space, so that it 
could be used for space navigation and positioning in the Solar System and 
beyond. In practice the initial event of the self-positioning process is used as 
the origin of the reference, and the axes are oriented according to the positions 
of the distant sources; all subsequent positions will be given in that frame. If 
one wants to further position the whole section of worldline of the receiver in 
some other external reference framj^ the location of the initial event in the 
external frame has to be known by other means. We describe our approach 
in flat Minkowski spacetime, and discuss the validity of this and other ap- 
proximations we have considered, for actual physical situations. The paper 
is organized as follows: after resuming the basic principles of null frames in 
Section [2l in Section [3] we show how localization takes place in our approach; 
then, in Section H] we discuss the sources of error, |5] is devoted to the discussion 
of some practical problems one has to face, while Conclusions are drawn in 
Section [Hi 



2 The Basic Null Frame and Definition of the Grid 

Let us consider a number of sources of periodic electromagnetic signals, at 
rest at spatial infinity, in a four-dimensional Minkowski spacetime. For our 
purposes, at least four sources are needed. Each of these sources is character- 
ized by the frequency of its periodic signals and by their directions in space; 
since the sources are supposed to be far away (i.e. at spatial infinity), their sig- 
nals can be seen as corresponding to plane waves. In the inertial frame where 
the sources are at rest, once Cartesian coordinates are chosen, we associate to 
each source a null four-vectoiO] / whose Cartesian contravariant components 
are given by 



T being the (proper) signal period, and n the unit vector describing the di- 
rection of propagation in the given frame. If in the same reference frame we 
consider the position four-vector 



For instance, in the Solar System events can be described by making use of the 
Barycentric Celestial Reference System (BCRS). 

^ Arrowed bold face letters like x refer to spatial vectors, while boldface letters like 
/ refer to four-vectors; Greek indices refer to spacetime components, while Latin 
letters label the sources. 






(2) 
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with respect to an arbitrary and yet unspecified origin, then we can define the 
scalar function X at the spacetime event identified by the position four-vector 
r 

Xir) = f r, (3) 

where dot stands for Minkowski scalar product. The scalar X might be thought 
of as the phase difference of the wave described by / with respect to its value at 
the origin of the coordinates. Four linearly independent four-vectors constitute 
a basis, or a fra me: we may think of choo sing four null four- vectors to serve as 



a basis (see e.g. iBlagojevic et al.l ( 120021 )). so that the four wave four- vectors 



{f{a), f(b), f{c), f{d)} in the form ([T]) constitute our null frame, or null tetrad. 
Then, accordi ng to the g enera l approach to coordinate systems and frames 



developed by ICoU et. all ( 120091 ) in connection with positioning systems, the 



four phase differences 

X(N) = f{N) r, N = a,b,c,d (4) 

defined at any event r whose coordinates are defined by with a,b,c,d 
labeling the sources, are null coordinates: in other words, they are spacetime 
functions with null spacetime gradient and, hence, they define a null coordi- 
nate system. Furthermore, the |X(7V)| are emission coordinates, since they are 
physically related to the reception of electromagnetic signals emitted by the 
sources. 



In detail, according to the general tetrad formalism (see e.g. IChandrasekhar 



( 1l983l )). the null tetrad allows the definition of the symmetric matrix 

'n{M){N) = f(M) ■ f{N), (5) 

which, in this case, has constant components, and whose inverse is determined 
by the relation 

^(M)(P)r/(^)(^) = 6[Z]. (6) 
Tetrad indices N = a,b,c,d are lowered and raised by means of the matrices 
V{M){P) and T]^^^^^^\ We can write the position four-vector r in the form 

r = xW/(;v) = X(^)/(^\ (7) 

and, as a consequence we see that the phase differences -^(at) are the compo- 
nents of the position four-vector with respect to the vectors 

fW = r^(^)(^nf^^^^^ (8) 
or, differently speaking, the functions 

XiN) = ^(^)(^)X(M) (9) 

are the components of the position four-vector with respect to the null tetrad 
vectors f{N)- It is useful to remark that while the frame {/(a), f{b), f{c), f{d)} is 
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constituted by null vectors, the frame {f^°'\f^''\f^'^\f^'^^} is constituted by 
space-like vectors. 



In summary, if we consider the hyperplanes conjugated to the n ull frame 
{.f(n.'), f(b), f{c), f{d)} vectors, we are able to define a spacetime grid (see lTartaglia 
( I2OO9I )). in which each event is identified by the relative phase of the electro- 
magnetic signals with respect to an arbitrary origin and, in this frame, the 
coordinates of each event are given by the functions Ix*^^-*!; equivalently, 

the phases |X(jv)| are the coordinates with respect to the space-like frame 



3 Localization within the Grid 



After having shown how to build a grid, we want to focus on how localiza- 
tion can be achieved within the grid. In particular, we suppose to deal with 
periodic signals in the forme of electromagnetic pulses, such as those coming 
from pulsars, and that these signals can be thought of as locally plane waves. 
Furthermore, we suppose that the user is equipped with a receiver able to 
recognize and count the pulses coming from the various sources, and a clock, 
that can be used to measure the proper time span between the arrivals from 
each source. 

Let us start with a toy model, where the emission from the sources is 
continuous and the phases of any signal can be determined with an arbitrary 
precision, at any event. We choose a starting event, from which the phases 
of each signal are measured, which is the origin of our coordinates (in other 
words, the event with r = 0, according to what we have described above), 
and three directions in space, defining the Cartesian axes of the inertial frame 
of the sources. We point out that even though the starting event is arbitrary, 
in order to correctly define the null frame, the directions of the sources in the 
sky have to be known: in other words, we have to know the unit vectors n 
for each source (and their proper frequency u too), which also enable us to 
calculate the matrices f]{N){M),V^'^^^^^'' of the given frame. 

To a subsequent event r, we associate the measured phases 

X(N) = /(TV) ■ r, (10) 

and, according to eq. (171), it is then possible to obtain the coordinates of the 
event r, in terms of the measured phases: 

r = X(„)/(") + X(,)/(^) + X(,)/('=) + X(,)/('^). (11) 

and to reconstruct the user's worldline. 
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Coming to a more realistic situation, such as the one in which the emitters 
are pulsars, we should consider that the received signals consist in a series of 
pulses and are not continuous. In this case, we may proceed as follows. First, 
we call "reception" the event corresponding to the arrival of a pulse from one 
of the sources. As a consequence, the position in space-time of an arbitrary 
reception event can be written in the form 

r = X(„)/(") + X(,)/W + X(,)/W + X(,)/(^), (12) 

with 



X(a)=n(a) +P, (13) 

X(b) = + g, (14) 
X(c) = n^c) + s, (15) 
X(^d)=n(d) + U!. (16) 

The X(7v)'s in the case of a continuous signal would be phases. Here they are 
given by an integer n(Ar), numbering the order of the successive pulses from 
a given source, and a fractional value: e.g. p means a fractional value of the 
cycle in X^"'\ and the equivalent holds for q, s, w, where < p,q, s, w; in eqs. 
f|T3l) - f|T6|) . only one of the p, q, s, w will in general be zero: when for instance a 
pulse from source "a" arrives p will be zero, while q, s, w will not; when a pulse 
from source "b" comes q will be zero and the other three fractions will not; and 
so on. Once we choose an arbitrary origin, we may count the pulses in order to 
measure the n(Ar), but we have no direct means to measure the fractional values 
p,q,s,w. However a procedure to determine these values can be obtained, 
based on geometric considerations: we suppose that the acceleration of the user 
is small during a limited series of reception events, so that we may identify 
the user's worldline with a straight line; furthermore, we also suppose that 
by means of his own clock the user can measure the proper time interval Tij 
between the i-th and j-th arrivals. With these assumptions we can proceed 
as follows to determine the fractional values p,q,s,w. Let us consider two 
sequence^ of arrival times from the sources; we have eight events, each of 
them in the form 



X(,),/('^) + X(,),/(^) + X(,),/('=) + X(,),-/('^), J = 1, .., 8, (17) 



where are expressions like flTBHT^ . The events are arranged in such a 

way that ri is the generic arrival of the signal from pulsar "a" , r2 is the arrival 
of the first signal of pulsar "b" after ri, is the arrival of the first signal of 
pulsar "c" after ri, and is the arrival of the first signal of pulsar "d" after ri 
(the pulsars are ordered from largest ("a") to shortest ("d") period); is the 



^ They may be subsequent or not, provided the total time span does not spoil the 
hypothesis of linearity of the worldline. 
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arrival of the next signal from pulsar "a" , and so on. The flatness hypothesis 
allows to write the displacement four-vector between two reception events in 
the form 



ri 



Indeed, the assumption that the worldline of the receiver is straight during a 
limited number of pitches of the signals can be used also to provide further 
information. In fact, let us consider three successive reception events i,j,k; we 
have 



The straight-line hypothesis allows us to write 



(19) 



AX, 



AX, 



(b)ji 



AX, 



AX, 



id)jt 



AX, 



{a)kj 



AX, 



ib)kj 



AX 



(c)kj 



(20) 



where Tji, t^j are the proper times elapsed between the i-th and j-th, and j-th 
and k-th reception events, respectively. These relations enable us to obtain 
the values we are interested in: in fact, we may arrange the coefficients of eqs. 
( IT71) in an 8 X 4 matrix (8 events, 4 sources): 



X, 



iN)i 



( nS") 


ni' + qi 


n\' + Si 


(d) , 

nl + wi 


ny +P2 


(b) 

ny 


ny + S2 


id) , 

ny +W2 


(1) , 
ny +P3 


(b) , 

ny +q3 


ny 


id) , 

ny +W3 


ny +pa 


(b) , 

ny + q4, 


ny + Si 


(d) 

nl 


(a) 


(b) , 

ny +q5 


ny + S5 


id) , 

ny +W5 


(i) , 
ny +pe 


(b) 

ny 


ny + se 


id) , 

ny +we 


7^7 + P7 


(b) , 

ny + qj 


n) 


id) , 

ny + wj 


\nl +ps 


(b) , 


ny + ss 


Ad) 

rig 



(21) 



As it can be seen, the p, q, s, w are zero along the main diagonals of the upper 
and lower square half matrices forming the whole matrix. This happens in 
correspondence of the arrivals of the pulses from the various sources: on the 
arrival of a pulse from "a" the corresponding p is zero, from "b" q is zero, 
and so on. Then, on using relations like ( |20|) we obtain the fractional values 
in terms of observed quantities, i.e. proper time intervals measured by the 
observer. For instance, we have 
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Pi 



= 0, (22) 



S3 



?4 



T62 

= - n^' - in? - nf) ^, (23) 

V / 7-73 



-n^f^-nf^-inf^-nfY-^, (24) 

4'^)_4'^) + (4'^)_4'^))Z2i^ ^^ = 0, (25) 
^ ^ T51 

^r-nr + ^^'-n?)^, (26) 



T"73 

nr-n?) + (nr^-ni^))^, (27) 

+ (4«)_^(«))!3i^ (28) 
^ ^ T51 

762 

0, ^3 = - nf + (n?^ - nf ) ^, (30) 

^ ^ T84 

P4 = nS«)-4«) + (nr-nn^, (31) 
^ ^ T51 

n<"-nf -(„<"- 4")^, (32) 

S4 = n^^ - ) - (n? - 4^) ^, «;4 = 0, (33) 
^ ^ 7-73 

and so on. Moving the pair of sequences and repeating the elaboration step by 
step, we are able to reconstruct the whole worldline of the receiver, in terms 
of measured quantities, i.e. proper times. 



S2 

P3 = ' - ^3 

<ls-n?-nr-{ni^-4^)^, (29) 



4 Biases and uncertainties 



In this Section we will discuss the sources of error and their importance 
in the positioning process that we have outlined. Before going into details of 
our analysis, we would like to preliminarily discuss the nature of errors and 
inaccuracies that affect our positioning process, in order to clarify how they 
can be dealt with to improve accuracy. 

Roughly speaking we can distinguish between systematic errors and un- 
certainties (or fluctuations). Systematic errors originate from mismodelling 
of the physical processes we are dealing with or from poor knowledge of the 
system parameters; they globally affect the positioning process and are ei- 
ther constant or time dependent. For instance, errors in the angular positions 
of the sources in the sky as well as in the periods of the pulses are system- 
atic errors; furthermore, if we deal with pulsars, their angular positions may 
change because of their proper motion and their periods because of energy 
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loss. As for the uncertainties, they are related to the stochastic variations of 
the quantities involved in the positioning process: for instance in the proce- 
dure of measuring the arrival times of the pulses, fluctuations are due to the 
detection device and to the user's clock as well as to the emission mechanism 
at the surface of the star; in principle also the turbulence of the interstellar 
plasma could have a role. Systematic errors produce global consequences with 
an unknown distortion of the re-built spacetime trajectory of the user; the un- 
certainties transform the worldline in an uncertainty stripe across space-time. 
Systematics can be reduced by improving our model and the knowledge of its 
parameters; a statistical analysis and the best technologies can help to reduce 
the impact of random disturbances. 

4.1 Model Limitations 

Let us begin by discussing the validity of our model, which is based on the 
propagation of electromagnetic signals originating from pulsating sources at 
rest in a given reference frame, in Minkowski spacetime. The question we would 
like to address is: how realistic is this model? Can it be used for positioning 
in actual physical situations? We start with some general observations. 

We explicitly work in fiat spacetime, thus eliminating the effects of the 
gravitational field. It is however obvious that, for positioning in the Solar 
System, the gravitational field of the Sun (and of the other major bodies) in- 
fluences the propagation of electromagnetic signalsQ The dimensionless mag- 
nitude of the static gravitational field of the Sun is of the order of 6q ~ 
'^J^® ~ 10~^ ^'^' ^IZI, and reaches its maximum value near the Sun, where 
Sq ~ 10^^. This field produces effects on the times of arrival of the pulses 
which are relevant for our purposes only if they change in a time comparable 
with the integration times used for our algorithm; only the radial component 
(with respect to the Sun) is important. The effects due to the motion of a 
user who travels with speed v in the radial direction r (and thus experiences 
a time varying gravitational field) in a time span 6t are expressed in terms of 
apparent fractional change of the period of the sources as Sq^^ ~ ^^^v • r 

6 X 10^*^ (^ ^ ( 30 km/s ) y"^*^© This fixes an upper limit to the acceptable 
6t in order the effect to be compatible with the required tolerance. One more 



^ En passant, we notice that since in our procedure only ratios of proper times 
measured in practice at the same spatial position are involved, the effect of the 
gravitational field on the user's clock can be safely neglected. 

^ The distance d is here expressed in astronomic units (A.U.) so that a distance 
equal to 1 corresponds to the position of the Earth. 

^ The velocity is in units of the average speed of the Earth along its orbit (30 km/s) 
and the rate of change is per year (y~^). 
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problem is that when the hne of sight to a source grazes the Sun the time of 
flight of the signals depends on the geometric curvature of the rays and on the 
Shapiro ti me delay, wh i ch is indeed huge depending on the apparent impact 



parameter (jStraumanru . |2004| ). However these disturbances can be dealt with 
either choosing sources which are located in the sky so that their lines of sight 
are far away from the gravitating bodies; or having recourse, which is appro- 
priate for many other reasons also, to a redundant number of sources, so that 
only the best are used each time. 

Another systematic error comes from considering the sources as being at 
rest, at spatial infinity, in a given reference frame, which is a very idealized 
situation. In fact, actual sources, such as galactic pulsars, have both a proper 
motion and a finite distance from the observer. Taking in to account estimates 



of the proper motion of real pulsars (IHobbs et al.l . l2005l ) it is possible to see 



that the rate of change of the angular positioij^is of the order 10^^ ( ^^^d ) 
per year. In practice these figures tell us that we are allowed to keep the 
position nominally fixed for months before correcting the value of the direction 
cosines, which is possible because the behaviour of the sources is known. As for 
the radial proper motion of the pulsar, which is commonly far worse known 
than its transverse motion, it is indeed contained in the times of arrival of 
the pulses, but can reasonably be considered as fixed during very long times, 
practically not affecting the positioning process. Finally it is known that the 
periods of the pulses are generally increasing with time, due to rotation energy 
loss of the star. However the decay rate is in general known, so that we may 
deal with this kind of change in the same way as we do for the position in the 
sky: keep the period fixed for a time compatible with the required accuracy 
(which would in any case be rather long), then introduce a corrected new 
value (which is known once the pulsar is given). As for glitches (sudden and 
random jumps in the frequency) they can be cured thanks to the redundancy 
of the number of sources above the minimum of four, corresponding to the 
dimensions of spacetime. 

In the case of pulsars there are many other effects that co uld affect the 



arriva l times of the signal; a thorough analysis can be found in iKramer et al. 



(120041 ). Besides what we have already mentioned, we should, for instance, take 
into account the emission mechanism at the surface of the star and the propa- 
gation across the interstellar plasma: these effects are responsible for the fact 
that each single incoming pulse is in general different from any other. The 
method to deal with this is the same adopted at radiotelescopes: the acquisi- 
tion of the data runs for a high number of pulses, then the signal is analyzed 
by folding, in order to extract the fiducial sequence with the desired accuracy. 



Actually, 


proper 


Straumann f 


20041. 



time variation of the period is of the order 10 I ^'^^^'^ ) y ^6t. 
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In the case of the positioning far shorter integration times are needed, because 
the pulsar is already known and the aim is simply to recognize it, rather than 
investigate on it. For completeness we should finally mention the problems 
related to the acquisition and elaboration chain in the receivers, however to 
discuss the technological aspects of the antenna and acquisition apparatus is 
beyond the scope of the present work. 

Finally, we would like to focus on the approximation that is implicit in 
our method for the conversion of a sequence of times of arrival of discrete 
pulses into the coordinates of the receiver. We exposed the method in Section 
131 where we limited the extension of the time series of the signals we used 
at each step of the process to an interval allowing to locally approximate 
the user's worldline with a straight line: how far is this hypothesis tenable? 
Given the user's clock accuracy St 0, we can define the maximum proper time 
interval Armax that can be considered in order to be self consistent with the 
straight-line hypothesis. Developing the worldline of the receiver in powers of 
its proper time up to the second order we see that, if a is the order of magnitude 
of the user's acceleration, and v his velocity, the following condition should be 
satisfied: 



For instance, if the user is moving in flat spacetime with 6t ~ 10^^° s, 
f = 5 X 10^ m/s and an acceleration a = 1 m/s^, we have Armax = 10~^ 
s, which corresponds to several periods of millisecond pulsars: enough both 
for the averaging away of the fluctuations of the single pulses, and for the 
piecewise reconstruction of the worldline. Actually, the deviation from the lin- 
earity of the user's worldline can also be due to the curvature of spacetime, 
i.e. to the presence of the gravitational field. We can give a similar estimate of 
the corresponding maximum proper time interval Armax by setting a = |V$| 
in flMl) where $ is the gravitational potential. For instance, for a = 10~^ m/s^, 
which is the order of magnitude of the gravitational field of the Sun at 1 A.U., 
we get for v = 10^ m/s, Axmax — 10~^ s. So, we see that there are actual 
physical situations where the hypothesis of linearity holds, and the procedure 
that we have described is meaningful. 



^ Into this 6t we should actually include also the drifts due to the proper motion 
and the period decay of the pulsar, which we keep constant during one step of the 
process; the latter are however expectedly far smaller effects than those due to the 
acceleration of the receiver. 




(34) 
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4-2 Errors in the procedure of position determination 



After having discussed the systematic errors related to the physical model 
underlying our procedure of position determination, we would like to turn to 
the analysis of the errors in the procedure itself which, as we are going to see, 
are connected both with our knowledge of the system parameters and with 
the measurement process. 

Actually, the coordinates of an event r are determined by solving eq. ffTOj) . 
We notice that eq. ffTOj) can be written in the form 

Ax = y (35) 

wher^ 





( 

/(a) 


-f{a) 


f2 f3 \ 

J (a) J (a) 




( 






A = 


J{h) 
J(c) 


-fib) 
-f{c) 


-f{b) -f{b) 
-f{c) -f{c) 


, X = 


X^ 

x^ 


> y = 


X(b) 

X(c) 




\f{d) 


-fid) 


-f{d) -ffd)) 




^X^ J 




\^{d)) 



with x° = ct, x^ = X, x"^ = y, x^ = z. Eq. ( l35l) is a linear system where the 
unknown vector x is obtained in terms of the matrix A (which has to be non 
singular) and the vector y: 



X = A-^y. (37) 

The entries of matrix A are related to the signal periods T(7v) = (c/[Jv))^^ 
and the direction cosines n^^-, = cT(^N')f(^j^--^ defining the angular positions of 
the sources; the fact that the n(Ar)'s, = a,b,c,d are different ensures that 
the matrix is not singular. As for the vector y, its entries are the measured 
phases determined (at each event) by the procedure described in Section |3l 

The entries of the matrix A are affected by systematic errors, while phase 
measurements have random errors. It is then possible to evaluate the maximum 
relative error in the solutions x of the linear system fl37j) by estimating the 
maximum errors on the entries of A, while a covariance analysis allows to 
quantify how random errors in phase measurements translate into errors in 
the determination of the coordinates of the spacetime event r, on the basis of 
the geometric properties of matrix A. 

^ We use X to refer to generic vectors belonging to M^, not to be confused with the 
four-vectors. 
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I n fact, after defining a suit able nor the following relation holds (see 
e.g. iDemidovic fc MaronI ( 1l966l )) 



6x < k{Af6A 



(3J 



where 6x 



Ml ' 



||AA|| 

ll^ll 



and k{A) = \\A\\\\A ^|| is the condition number 
of the system ( l35l) . If we suppose that the relative errors in periods AT/T, 
and direction cosines An* are roughly the same for all sources in all directions, 
then eq. (138!) allows to conservatively estimate the relative error in the form 




6x < k{Af 

We remember that it is always k{A) > 1; in our case, we may write 

1 ^ 1 

E 



(39) 



k{A) oc 



|det(A)|^ 



rp2 ^2 



N=a,b,c,d (N) 



(40) 



As a consequence, we see that the relative error is minimized when: (i) the 
determinant is maximized and (ii) periods are minimized; since the spatial 
components of A are the direction cosines, the determinant is maximized 
when the volume spanned by these directions is maximized. To fix ideas, for 
AT/T ~ 10-^ An ~ 10"^ it is 5x < k{Af 10-\ 



Further insight on the accuracy achievable by our positioning procedure 
can be obtained by discussing the geometric properties of the system (1371) . 
which ultimately depends on the actual angular positions of the sources in the 
sky: we refer to the Geometric Dilution Of Precision (GDOP), which is based 
on the covariance matrix of the errors in position determination, and provides 
a measure of how fit the set of sources is: indeed, it is related to the geometric 
properties of the matrix A and a great accuracy can be obtained if sources 
are chosen that are sufficiently scattered i n the sky (this is pretty rauch li ke 
what happens with GPS satellites, see e.g. iHofmann-Wellenhof et al.l ( 120011 )). 
The covariance matrix of the positioning errors, as determined by the phase 
measurements, can be expressed as 



cov X 



A^ (cov y) A 



-1 



(41) 



If the phase measurements errors from each source are uncorrelated and Gaus- 
sian distributed, with z 
matrix turns out to be 



sian distributed, with zero mean and variances cr^^^, the position covariance 



cov X 



A diag((7(„),(T(';,),(T(-^), 



cr, 



-1 



(42) 



10 



For instance: 11^ I 
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The GDOP is defined by 



GDOP = ^Tr(cov x) (43) 

To understand the meaning of GDOP, let us consider the case when all phase 
measurements have the same variance cr^, then cov y = diag(cr^, cr^, cr^, o"^); 
furthermore it is {A^A)~^ = ^^J^^^2 (cofA)'^ (cofA) so that we can write 



GDOP = v'Tr(cofA)^(cofA) (44) 

Again, we see the role of |det(A)| in minimizing the errors: GDOP is minimized 
when the solid angle spanned by the directions of the sources is maximized. 

The error in position determination can be minimized by using a number 
n of sources greater than four, which will result in an overdetermined linear 
system in the form fl35|) , Bx = y, provided that now is an n x 4 matrix and 
y is a vector of n measured phases. Redundant equations allow to improve 
the accuracy of t he so lutions via least squares estimation techniques (see e.g. 
Stark fc Wood] fll986h ). so that Eq. (gl]) is generalized tcEi] 



cov X 



{B^B)-^B^ (cov y) B{B^By^ (45) 



If, as before, we consider the phase measurements errors to be identically Gaus- 
sian distributed with the same variance and independent, we obtain equation 
pri) back, also in the case n > 4: in summary, we may say that the components 
of {B^B)~^ determine how phase errors translate into errors of the computed 
spacetime position. 



5 Technical problems 



The purpose of the present paper is not to design a working navigation 
system, rather to describe a new method allowing relativistic positioning when 
a set of pulsed sources of electromagnetic signals and an onboard clock are 
available. We refer mostly to pulsars because they represent natural sources 
of regular pulses and they are suggestive of the old navigation at sea using 
the stars; an advantage of using pulsars is the fact that they can practically 
be considered as being at infinite distance and occupying fixed positions in 



the sky I. however a major drawback is the weakness of their signals. Being 



aware of this, we have gone on with our example of pulsars without facing the 



Notice that if n = 4 it reduces to Eq. ([^ . 
^■^ Actually this is not the case, but we have already commented on how to treat the 
issue in Section 14.11 
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practical issue of the detection of the signals, but our method can perfectly 
work with any other artificial source of pulses provided one knows the law 
of motion of the source in a given reference frame. For the use in the Solar 
system, one could for instance think to lay down regular pulse emitters on the 
surface of some celestial bodies: let us say the Earth, the Moon, Mars etc. The 
behaviour of the most relevant bodies is indeed pretty well known, so that we 
have at hands the time dependence of the direction cosines of the pulses: this 
is enough to apply the method and algorithm we have described and the final 
issue in this case would be the position within the Solar system. In principle 
the same can be done in the terrestrial environment: here the sources of pulses 
would be onboard satellites, just as it happens for GPS, but without the need 
of continuous intervention from the ground: again the key point is a very good 
knowledge of the motion of the sources in the reference frame one wants to 
use. 

In the previous sections we have many times mentioned the fact that the 
number of independent sources to be used is at least four. The reason for 
that number is that the space-time coordinates needed to localize the user 
are of course four (three space components plus time). There are however 
many reasons for using a bigger number of emitters. If we have n > A pulsars 
(or equivalent sources) we may apply our method to all possible quadruples 
contained in n: the position will be determined as an average of the results 
obtained from all quadruples. This would be one more way to average the 
effect of random disturbances at the emission of the pulses out. Furthermore 
if one of the sources fails or disappears for any reason (e.g. because it is 
eclipsed) the localization process is not interrupted; if one new source comes 
in, provided its position in the sky and proper recursion time are known, the 
sequence of its pulses is hooked to the main sequence of the arrival times and 
may be used for further positioning. What matters is that all arrival times 
from all visible sources are arranged in one single sequence identified by the 
onboard clock and that the main sequence is not interrupted. To say better, 
the maximum tolerable total blackout of the sequence should not last more 
than what can be reconstructed by extrapolation from the last portion of the 
world-line of the observer. With reasonable accelerations at the typical speeds 
of an interplanetary travel this may be as long as a few seconds. 

Of course everything works if the clock used by the observer is a good 
and reliable clock. We are considering atomic clocks fit for being carried on 
board a spacecraft. The problems with the clocks arise from both random 
and systematic instabilities. The former contribute to the general uncertainty 
in the positioning, whilst the latter introduce an apparent shortening of the 
period of the received pulses. A strategy to minimize the effects of the drift 
in the frequency of the clock is similar to the one adopted for GPS. It is 
redundancy: more than one atomic clocks, and of different types, should be 
carried by the observer (GPS sateUites have four clocks on board): the arrival 
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times would be given by an average of the different readings. Tlie redundancy 
would partly reduce also the effect of the drift of the frequency of the clocks, 
but this would emerge in long time. The accur acy of portab le present day 
atomic clocks is in the order of 1 x 10~ (see e.g. lAshbvl fl2003h ) which means 



that the measured repetition time of a millisecond pulsar would be significantly 
affected by the frequency drift of the clock over times of the order of months, if 
not years. It would in any case be advisable to periodically check the accuracy 
of the measurement from some reference station on the Earth, and the best 
strategy would be to check the position and compare with the one obtained 
by self positioning, rather than to try to directly verify the clock, which would 
be possible but it is a very delicate task indeed and presupposes, any way, 
a good knowledge of the relative motion between the space observer and the 
ground station. A periodic check up (after weeks or months intervals) would 
in any case be cheaper than a continuous guidance from the ground. 

Coming back to the initial lines of the present section, we stress again the 
fact that ours is not the design of a working device, but rather the description 
of a principle method to attain self positioning by means of the signals emitted 
by periodic sources from known positions in the sky. In the same spirit we 
have here shown that we are considering also specific strategies to tackle some 
relevant practical problems one would have to face when coming to implement 
our proposal in the real world. 



6 Conclusions 



We have described an operational approach to the use of pulsating signals 
for positioning purposes; in particular, pulsars signals can be used. Our pro- 
cedure is based on the definition of a null frame, by means of the four-vectors 
associated to the signals in the inertial reference frame where the sources are 
at rest (so that the emission directions and the frequencies of the pulsating 
signals have to be known) and far away (so that their signals can be dealt 
with as plane waves). The procedure is fully relativistic and allows position 
determination with respect to an arbitrary event in flat spacetime. Once a null 
frame has been defined, it turns out that the phases of the electromagnetic 
signals can be used to label an arbitrary event in spacetime. If the sources 
emit continuously and the phases can be determined with arbitrary precision 
at any event, it is straightforward to obtain the coordinates of the user and 
his worldline. However, actual sources emit signals which consist of a series of 
pulses: so, we have developed a simple method that can be used to determine 
the user's worldline by assuming that the worldline is a straight line during 
a proper time interval corresponding to the reception of a limited number of 
pulses, which means that the effects of the acceleration are negligibly small. 
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In discussing the source of errors which affect the positioning process, we 
considered systematic errors, due to mismodelhng or to the poor knowledge 
of the system parameters, and uncertainties, due to the stochastic changes of 
the quantities involved in the procedure. While the former introduce global 
consequences resulting in an unknown distortion of the reconstruction of the 
user's worldline and can be reduced by improving our model and the knowledge 
of its parameters, the latter transform the worldline in an uncertainty stripe 
across spacetime and can be reduced by means of a statistical analysis and 
an improvement of the technologies involved. We have discussed the effects of 
the gravitational field, of the motion of the sources and of the variability of 
their emission, the validity of the hypothesis of linearity, and we have shown 
that their impact is small enough or can be corrected for a number of actual 
physical situations, so that the procedure that we have described can be safely 
applied. Also, we have evaluated the accuracy that can be achieved, by taking 
into account the maximum errors in the parameters of the system and by 
quantifying the impact of random measurement errors on the determination 
of the spacetime trajectory, introducing the Geometric Dilution Of Precision. 

We did not enter into several technological issues that are important for 
building an actual positioning system based on our procedure, for instance the 
extreme weakness of signals from pulsars and the required design and sensitiv- 
ity of the detectors. These and other issues require further investigations and 
perhaps technological improvements, however we believe that this approach 
could be useful for defining an autonomous and relativistic positioning sys- 
tem in space, which ultimately transfers the basic positioning frame from the 
Earth to spacetime, according to a truly relativistic viewpoint. 
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